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Following the TMD evolution scheme recently proposed for the unpolarized and the Sivers 
distribution function, we propose a simple strategy to take into account this TMD Q 2 
dependence in our phenomenological extraction of the Sivers function from SIDIS data. 
New results are presented and possible future applications are discussed. 

The exploration of the 3-dimensional structure of the nucleon, both in momentum and in 
configuration space, is one of the major issues in high energy hadron physics. Information on 
this 3-dimensional structure is embedded in the Transverse Momentum Dependent distribution 
and fragmentation functions (TMDs). The Sivers function, which describes the number density 
of unpolarized quarks inside a transversely polarized proton, is particularly interesting, as it 
might provide information on the partonic orbital angular momentum. 

So far, all phenomenological fits have either neglected the QCD scale dependence of TMDs 
(which was unknown) or limited it to the collinear part of the unpolarized PDFs, according 
to the DGLAP evolution. Here, we present the first attempt to take into account the TMD 
evolution as proposed by Aybat, Collins, Qiu and Rogers [TJ[51|3] in the analysis of the Sivers 
asymmetry data and show how these new results compare with the previous extractions. Even- 
tually, such a scheme will provide a complete TMD factorization framework for a consistent 
treatment of all SIDIS data. 

In Ref. [3] we showed how the QCD evolution equation of the TMDs in the coordinate space 
proposed in Refs. [2] and [3] can be expressed in a simplified way, taking the renormalization 
scale /i 2 and the regulating parameters £p and £d all equal to Q 2 , as 

F{x > b T ;Q) = F(x,b T ;Q )R(Q,Q ,b T ) exp |- 5jf (6 T ) In ^| , (1) 

where F can be either the unpolarized parton distribution, F(x,bx',Q) — fq/ p (x,bT',Q), the 
unpolarized fragmentation function F(x,bT',Q) = Dh/q(z, for; Q), or the first derivative, with 
respect to the parton impact parameter 6 T , of the Sivers function, F(x, b T \ Q) = f' 1T *(x, b T ; Q)\ 
9kQ>t) is an unknown, but universal and scale independent, input function, while R(Q, Qo, br) 
is the evolution kernel 

R(Q, Qo, br) = cxp {in £ T ^ 7 aV) + P ~1f (h, %) ) • (2) 

Qo Jqo \J J Qo m V ^ J \ 
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The anomalous dimensions jp and jk appearing in Eq. |2]), are given, at order 0(a s ), by [2] 



7f(m; ^o) = ai a (ji) — f ~ -In—") 7k(m) = 2 ^ ' ( 3 ) 

/l 2 7T \2 jLT 



The Q 2 evolution is therefore driven by the functions gxibr) and R(Q,Qo,br)- While the 
latter, Eq. can be easily evaluated, numerically or even analytically, the former, is essentially 
unknown and will need to be taken from independent experimental inputs. 

The appropriate Fourier transforms allow us to obtain the distribution and fragmentation 
functions in the momentum space: 

1 f°° 

f q /p(x, kj_; Q) = — db T b T J (kj_b T ) f q / p (x, b T ; Q) , (4) 



2tt 

D h / q (z,p x ;Q) = ^- / d6 T b T J (k T b T ) D h/q (z,b T ;Q) , (5) 
Z7r Jo 



— 1 f°° ~ 
f x 4 {x,k±;Q) = / d&r & r Ji(fc±6r) f 1T q {x, b T ;Q) , (6) 

27r/e^ Jo 

where Jo and J\ are Bessel functions, while f q / p is the unpolarized TMD distribution function 
for a parton of flavor q inside a proton, D h / q is the unpolarized TMD fragmentation function 

for hadron h inside a parton q and /{j? is the Sivers distribution describing unpolarized partons 
inside a transversely polarized proton. 

The unknown input functions <7k(&t) and F{x 1 bx' 1 Qo) inside Eq. |lj have to be appropri- 
ately parameterized. As already anticipated, gxibr) is a non-perturbative, but universal func- 
tion, which in the literature is usually parameterized in a quadratic form: <7^(6t) = \ gi b\. As 
in Ref . [3J we will adopt the results provided by a recent fit of Drell-Yan data [5] , and assume 
<?2 = 0.68 GeV 2 . The input functions F(x,bx',Qo) are parameterized by requiring that their 
Fourier-transforms, which give the corresponding TMD functions in the transverse momen- 
tum space, coincide with the previously adopted fc^-Gaussian forms, with the x dependence 
factorized out. As shown in Ref. [3], one finds 

J q/P (x, b T ; Q) = f q/p (x, Qo) R(Q, Qo, b T ) exp j-6 2 , (a 2 + ^ In |Q j (7) 
D h/q (z, b T ; Q) = ^D h/q (z, Qo) R{Q, Qo, br) exp |-&| (V + | In ^ J (8) 
f£(x,b T ;Q) =-2 7 2 / 1 L T (x;Qo)i?(Q,Qo,6T)6Texp|-6^ 7 2 +f ln ^) } ( 9 ) 

with a 2 = (fci)/4, /3 2 = (pi>/(4z 2 ), 4 7 2 = = and in Eq. (g. 

Eqs. (J^J) - (JsJ) show that the Q 2 evolution is controlled by the logarithmic Q dependence of 
the bx Gaussian width, together with the factor R(Q, Qo, br)- for increasing values of Q 2 , they 
are responsible for the typical broadening effect already observed in Refs. [2] and [3]. 

As R(Q,Qo,bT) shows a weak dependence on (large) bx (i-e. small fcj_), we can assume 
R(Q,QQ,br) to be constant in br and compute the Fourier transforms of the evolution equa- 
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Table 1: \ 2 contributions corresponding to our three fits, for some experimental data sets of 
HERMES and COMPASS experiments. 





TMD Evolution (exact) 


TMD Evolution (analyt.) 


DGLAP Evolution 




Xlt = 255.8 


Xlt = 275.7 


Xtot = 31 5 - 6 




Xlo.f = 1-02 


xL./ = i-io 


Xio./ = 1-26 




Xl = 10.7 


Xl = 12.9 


Xl =27.5 


HERMES 7T+ 


Xl = 4.3 


= 4.3 


X 2 =8.6 




X?V = 9.1 


X 2 Pt = 10.5 


X 2 Pt = 22.5 




Xl = 6.7 


Xl = H.2 


= 29.2 


COMPASS h+ 


Xl = 17.8 


X 2 = 18.5 


xl = 16-6 




X 2 Pt = 12-4 


X 2 Pt = 24.2 


x% T = n.8 



tions 0, ([8]) and ^ analytically, to find 

e -fci/- 2 

f q /p(x, k±;Q) = f q /p(x, Q ) R(Q, Q ) (fO) 

2 / 2 

D h / q {z,p±; Q) = D h/ (z, Qo) R(Q, Q ) g~ ( 1]L ) 

' ' TTW F 

- fci (fc 2 ^ 2 e - k l/ w s 

& N f q /pt {x, k ± ;Q) = jj- i V2~e A N f q/p r(x, Q ) R(Q, Q ) ^ , (12) 

where f q / p (x,Qo) and D h / q (z,Qo) are the usual integrated PDF evaluated at the initial scale 
Qo, and A N f q / p t (x, Qo) gives the x dependence of the Sivers function [1] [—(2k±/M p ) f^ T = 
A. N fq/pt]- Most importantly, w 2 , Wp and w% are the "evolving" Gaussian widths, defined as: 



w 



■2 _ -1.2 \ o „ 1„ Q ...2 _ /„2 \ , o_2„ i_ Q ...2 /i.2\ , o- i_ Q 



= (fci) + 2 52 In ^ , w F = (pi) + 2z 2 g 2 In ^ , w 2 s = (fcjj s + 2g 2 In ^ • (13) 



Qo 



Notice that the Q 2 evolution of the TMD PDFs is now determined by the overall factor R(Q, Qo) 
and, most crucially, by the Q 2 dependent Gaussian width w(Q,Qq). 

It is interesting to point that the evolution factor R(Q, Qo), controlling the TMD evolution, 
is the same for all functions (TMD PDFs, TMD FFs and Sivers ) and is flavor independent: 
consequently it will appear, squared, in both numerator and denominator of the Sivers azimuthal 
asymmetry and, approximately, cancel out. Therefore, we can safely conclude that most of the 
TMD evolution of azimuthal asymmetries is controlled by the logarithmic Q dependence of the 
k± Gaussian widths w 2 (Q,Qo), Eq. (13 1. 

The aim of our paper is to analyze the available polarized SIDIS data from the HERMES 
and COMPASS collaborations in order to understand whether or not they show signs of the 
TMD evolution proposed in Ref. [3J. In particular we perform three different data fits of the 
SIDIS Sivers single spin asymmetry A™ (0h ~ 0s) measured by HERMES and COMPASS: a fit 
(TMD-fit) in which we adopt the TMD evolution equations of Eqs. 0-(|9]);a second fit (TMD- 
analytical-fit) in which we apply the same TMD evolution, but using the analytical approxima- 
tion of Eqs. ( fTo| ), ( 11 ) and ( 12 ); a fit (DGLAP-fit) in which we follow our previous work, as done 
so far in Ref. [HI [7], using the DGLAP evolution equation only in the collinear part of the TMDs. 
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Table I shows the main results of our fit- 
ting procedure. The best total x\of which 
amounts to 256, is obtained by using the 
TMD evolution, followed by a slightly higher 
Xtot °f the analytical approximation, and a 
definitely larger Xtot — 316 corresponding to 
the DGLAP fit. The difference of about 60 
X 2 -points between the TMD and the DGLAP 
fits is heavily concentrated in the asymme- 
try for 7r + production at HERMES and for 
h + production at COMPASS, especially when 
this asymmetry is observed as a function of 
the x-variable. It is important to stress that, 
as x is directly proportional to Q 2 through the 
kinematical relation Q 2 = xy s, the x behav- 
ior of the asymmetries is intimately connected 
to their Q 2 evolution. While the HERMES 
experimental bins cover a very modest range 
of Q 2 values, from 1.3 GeV 2 to 6.2 GeV 2 , 
COMPASS data raise to a maximum Q 2 of 
20.5 GeV 2 , enabling to test more severely the 
TMD Q 2 evolution in SIDIS. These aspects are illustrated in Fig. [I] where the SIDIS Sivers 
asymmetries A s ^^ h ~^ s ^ > obtained in the three fits are shown in the same plot. It is evident 
that the DGLAP evolution seems to be unable to describe the correct x trend, i.e. the right 
Q 2 behavior, while the TMD evolution (red solid line) follows much better the large Q 2 data 
points, corresponding to the last x-bins measured by COMPASS. 

In conclusions, we have analyzed the Sivers effect by up-grading old fits with the addition of 
the most recent HERMES and COMPASS SIDIS data, and by applying, for the first time, TMD 
evolution equations. We have compared the results obtained using TMD evolution equations 
with those found by considering only the DGLAP evolution of the collinear part of the TMDs. 
Our results give evidence that SIDIS data support the TMD evolution scheme, although further 
experimental data, covering a wider range of Q 2 values, are necessary to confirm this. 
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Figure 1: The results obtained from our fit of 
the SIDIS A s J^ Wh " 0s) Sivers asymmetries ap- 
plying TMD evolution (red, solid lines) are com- 
pared with the analogous results found by us- 
ing DGLAP evolution equations (blue, dashed 
lines). The green, dash-dotted lines correspond 
to the results obtained by using the approxi- 
mated analytical TMD evolution. The experi- 
mental data are from HERMES [8 (left panel) 
and COMPASS [H] (right panel) Collaborations. 
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